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Abstract. 

We present an analysis of the two-dimensional Schrodinger equation for two electrons 
interacting via Coulombic force and confined in an anisotropic harmonic potential. 
The separable case of uj y : uj x = 2 is studied particularly carefully. The closed-form 
expressions for bound-state energies and the corresponding eigenfunctions are found 
at some particular values of oj x . For highly-accurate determination of energy levels at 
other values of u) x , we apply an efficient scheme based on the Frobenius expansion. 
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1. Introduction 



The Hookean system composed of N interacting particles trapped in an external 
harmonic potential 

V(r) = ^[u 2 x x 2 + u 2 y 2 + u J 2 z 2 } (1) 

has been considered first as a model of nucleus in the early days of nuclear physics [I]. 
In the case the two particles interact via central forces, the Hookean system enjoys a 
particularly nice feature that the center-of-mass motion may be separated out. Because 
of its relatively easy tractability, the system has nowadays become a standard tool for 
testing quality of many-body approximation methods in studying correlations between 
electrons. Recently, the Hookean system attracted renowned interest, since the progress 
in semiconductor technology has allowed the formation of quantum dots (QDs), which 
are the systems of few electrons confined on the scale of hundreds of nanometers in an 
approximately harmonic potential [2]. For the simplest QD, consisting of two electrons 
(N = 2), the Hamiltonian can be written as 



i=i 1 

where m* is the effective electron mass and e* the effective dielectric constant. 
Introducing center of mass and relative coordinates 

R = ~(n + r 2 ) = (X, Y, Z), r = r 2 - n = {x, y, z), (3) 

the above Hamiltonian is separated into H = H + H r , where the center of mass 
Hamiltonian 

hR =^ + ™MX 2 + " 2 y Y 2 + uj 2 z Z 2 ) (4) 

is exactly solvable. The problem is thus reduced to the Schrodinger equation for the 
relative motion 

H r %l) = e r tp, (5) 
described by the Hamiltonian 

H r = -h^/m* + + uy + uy] + ( 6 ) 

where r = |r|. By a transformation x = x/a*,y = y/a*, the relative motion equation is 
rescaled to the form 

r d d d . o . o . o . o , o . o X 



dx 2 dy 2 dz 2 yy ^x 2 + y 2 + z 2 

where e r = R*e r , u x ^ tZ /2 = uj Xjy>z (h/ m^a 2 ) with the effective Bohr radius a* = 
h^Ex/m^e 2 and the effective Rydberg constant R* = m^e 4 /fi 2 e 2 . In the following, we 
skip the hats over spatial coordinates x,y,z, frequencies 0J xyz and relative energy e r . 
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In the case of some particular relation between the frequencies u x , oo y and u z , the 
problem of the relative motion becomes easier to solve, and at some specific values of 
frequencies it even appears quasi-solvable, i.e. one of the states of the spectrum can 
be analytically determined. The most studied example, important for describing QDs 
in a constant magnetic field, is the case the two of confinement frequencies are equal, 
e.g. uj x — uiy — lo. For such axially symmetric system there are three cases where the 
problem of the relative motion (J7j) is integrable: 1) uj z = u, being separable in spherical 
coordinates, 2) uo z : uo = 2, being separable in parabolic coordinates, and 3) u z : uo = 1/2, 
which is one of few examples of integrable problems that is generally not separable [3]. 
The quasi-solvability of the case 1) has been known for longtime [H [5], and recently it 
has been demonstrated [6] that also in two other cases, one of the bound-state solutions 
(not necessarily the ground-state) can be obtained in a closed-form at some particular 
values of uj. 

Another example, which is important in studying two-dimensional (2D) QDs, is 
that where one of the frequencies is much larger than the other two, e.g. uj z — > oo. 
The 2D approximation is usually justified in practice, since QDs are realized by lateral 
confinement of electron gas at the boundaries of semiconductor nanostructures. The 
above problem is integrable in two cases: that of uj x = ui y , and that of u y : u x — 2 
(which is equivalent to the case of u y : uo x = 1/2, under the exchange of uo x and u y ). 
In the first case, the harmonic potential has circular symmetry and the problem of the 
relative motion becomes separable in polar coordinates. This case was much studied 
in the literature [3, [8], [9] and the closed-form solutions for particular frequencies have 
been derived [10]. The influence of the anisotropy aspect ratio 0J x /uJ y on the QDs 
properties has been recently studied in Ref. pj], [12], where the separability of the case 
of uj y : uj x = 2 in parabolic coordinates has been pointed out. Here we shall study this 
integrable case in detail. The main result of our paper is the demonstration that for 
particular values of uj x the exact solutions can be found in a closed form. Although the 
property of quasi-solvability appears only at particular confinement frequencies, it is 
very appealing, as the knowledge of exact bound-state energies and the corresponding 
eigenfunctions in a closed form can be used for extracting information about the model 
and for testing the performance of approximation methods. Therefore, we discuss the 
quasi-solvable cases in detail, providing the values of confinement frequencies u x as 
well as the corresponding energy eigenvalues and the closed-form solutions for wave- 
functions. For determining the energy spectrum at arbitrary values of uj x , we develop 
a numerical scheme based on the Frobenius method (FM), which was successfully used 
for calculating bound state energies of one-dimensional [Hi HH [15] and spherically 
symmetric potentials [T6l [T7] . The comparison of the bound-state solutions determined 
within the FM with the exact ones allows us to demonstrate how good is the performance 
of the numerical approximation scheme. 

The outline of our work is as follows. In section [21 we give a brief discussion of 
the spectrum in 2D case of motion ([H]), paying a particular attention to the integrable 
case of ujy : uj x = 2. In section [3j the quasi-solvable cases are derived. In section [4[ the 
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calculation within the FM are described and a comparison with exact results is given. 
The paper ends with concluding remarks in section [5j 

2. Theoretical study 

In the case of laterally confined QD, the relative motion equation (JTj) takes the form 

+ u 2 x x 2 + u 2 y 2 + -]V(r) = £ r V(r), (8) 



L dx 2 Q y 2 V» p J 

where p = y^x 2 + ?/ 2 . The singlet (s) and triplet (t) spin states of the two-electron 
system correspond to symmetric and antisymmetric spatial wave-functions, respectively. 
As the centre-of-mass coordinate remains the same upon the interchange of electrons, 
the symmetry requirement reduces to the symmetry of the relative wave-function under 
inversion r — > — r . Because of the invariance of equation (jSJ) to reflections about the x- 
and y-axis, the (x, y)-parity of spatial wave-functions is well-defined. The parity (+, +) 
or(— , — ) corresponds to spin singlet eigenfunctions, and the parity (+, — ) or (— , +) to 
the spin triplet ones. 

The simplest case of circular symmetry has been much studied in the literature 
[3 13 |j. Here we discuss the other integrable case, u y : u x = 2. Introducing two- 
dimensional parabolic coordinates through 

x = vm,y= 2(^1 - vl),V2 > o (9) 

the Schrodinger equation is transformed into the form 

H r (vuV2)^(Vu V2) = £ V ip(Vu V2) (10) 
with the Hamiltonian given by 

ffr(, ' i ' ,fe)= ^rb h S^fi + r» ( '' ?+ '' 26) 

- \riWM + viM - + 2] (11) 

Now, it is easy to see that by setting u y = 2u x and ^(771,772) = git^i)^^) into 
(fi~T|) one gets two ordinary differential equations of the identical form 

~ £VT t + ^ + Pj]gj{j]j) = °' J = h 2 (12) 

that are coupled by the condition f3\ + $2 = 2 to be satisfied by the separation constants 
(3\ and Z^. The integr ability is explained by the fact that in the case of cu y : cu x = 2 an 
operator commuting with the Hamiltonian fill I) exists in the form [TT] 

A = {L z ,p x } + 2u 2 x yx 2 (13) 

P 

where the bracket {} denotes anticommutator. The operator A commutes with the 
x— reflection but anticommutes with the y— reflection operator, and its eigenvalue is 
given by 5 = (Pi — /3 2 )/2. Therefore, the eigenstates corresponding to 5 7^ 0, i.e. 
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f3\ 7^ fa, are doubly degenerate with respect to the sign of S, and those with 5 = 0, i.e. 
(3\ = (3%, are nondegenerate. 

The bound-states of the QD correspond to the values of e r and (3j such that the 
functions ^(771,772) vanish, as 1 771], 772 — > 00. The eigenfunctions of (fl2l) which vanish 
at infinity can be ordered by the quantum number rij that counts its nodal points on 
—00 < 77j < 00, and will be denoted by g^'irjj, j3j). The reflection symmetry of of ffTSl 
divides g( nj \r)j, fy) into two parity types, depending on rij being even or odd. As the 
inverse coordinate transformation 

rji = t^Vp + V and r/ 2 = y/p-y, 
\x\ 

shows discontinuity of in the 7/1 coordinate at positive y-axis, the condition emerges 
that the solution of fTTOT) be a product of functions of the same parity in 771 and 772, since 
only such products have the appropriate continuity properties in the whole space. In 
the case of 5 = 0, which corresponds to the functions g(r)j,/3j) having the same number 
of nodal points, the singlet state (s) with (x, y)— parity (+, +) is represented by the 
product of functions being even in 771 and 772, while the triplet state (t) with (— ,+) 
is represented by the product of odd functions. For 5^0 the solutions of ( flOl) with 
well-defined (x, y)-parity are easily constructed as 

^ ,n2 \vuV2) = g (ni) (m, Pi)g {n2) (m,P2) ±g (ni) (v2, fog™ fa, fa), (u) 

where the sign (+) corresponds to the singlet /triplet state with (x,y)— parity 
(+,+)/(—,+), and the sign (— ) to the triplet /singlet state with (+,—)/(—,—) in the 
case of n\ and ri2 being even/odd. 

In order to visualize the manifestation of degeneracies discussed above, we consider 
the behavior of bound-state energies for varying anisotropy uj x /uj y . The numerical 
calculation have been performed with the help of the linear Rayleigh-Ritz method using 
the basis of two-dimensional harmonic oscillator eigenstates ("exact diagonalization"). 
The energies of low-lying states of the relative motion, e v , are plotted in figured], for a 
fixed frequency uj x = 0.5 in function of the parameter u y . It can be observed how the 
states with spatial parity, (+, —),(—, +) become degenerate with those of (+, +), (— , — ), 
respectively, in the case discussed above, i.e. ui y : uj x = 2, and in the case equivalent to 
it after exchange x y, i.e. u x : u y = 2. The visible degeneracy of the spectrum at 
u y = uj x = 0.5 corresponds to a circular symmetry of the harmonic potential. 

In the following we will discuss the spectrum of the anisotropic QD in the case of 
uj y : u x = 2, which can be determined by solving a coupled pair of one-dimensional 
equations (fTSI) . A simple scheme based on the FM will be used for the numerical 
determination of the solutions in the above case. The exact solution will be provided 
for particular values of u x . 

3. Exact solutions 

Let us to come to main purpose of this paper which is the demonstration that the 
problem of a 2D QD, confined by the potential w x x 2 + uj^y 2 , is quasi-solvable for 




Figure 1. The behaviour of the relative motion energy levels of a QD at to x = 0.5 in 
function of ui y . The vertical lines correspond to the cases u> x : u> y = 2 and & y : u) x = 2. 
The circularly symmetric case of u/ y = u) x is marked by the broken line. 



uj y : uj x = 2, i.e. that some of the solutions of the relative motion equation (jSJ) can 
be found algebraically for certain values of uj x . We are looking for a solution of the 
form ^(771,772) = g(r]i, Pi)g(rj 2 , ^2), where Pi + P2 = 2, and g(r]j,Pj) fulfils the ordinary 
differential equation (fT2l . As the form of EqJT2]is similar to that of the sextic anharmonic 
oscillator, the quasi exact solutions can be represented as [TT 



^i,^)=e-^^apV^ (15) 

i=0 

where ag 7^ 0, and u = 0, 1 correspond to the even- and odd-parity case, respectively. 
Substituting (fI31) into (fl2|) and comparing the coefficients of like powers of r]j, one gets 
the three-term recursion relation for the expansion coefficients in the form 

Aa^, + B^ajp + Ciopi = 0, (16) 

where 



A = (2i + i/ + l)(2i + i/ + 2) 

Cj = e r + (1 - Ai - 2u)uj x (17) 

with af = for z < 0. It turns out that the QD equation (jHJ) is exactly solvable if the 
series in the representations (|15p for both g(t]i, Pi) and g(r)2, P2), terminate after a finite 
number of terms. This gives the three conditions 

£ r =0^(3 + 471 + 2^) (18) 
a^ +1 {uj x ,e r ,pi) = and a^ 1 (a; a .,e r ,&) = 0, (19) 
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where n denotes the number of the highest non-vanishing term which appears to be 
necessarily the same in both series. Using the recurrence relation ffl6|) . the coefficients 
of the series may be conveniently represented in the form 



•Pi 



1 



,n+l a 0') 



(2(n + 
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A x 
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A 2 







(2(n + 1) +u)\ 





















C n -i 


R (i) 


A n -i 














B U) 

J ->n 


J 



,3 = 1,2. (20) 



(i) 

By calculating the above determinant for a chosen value of n, the polynomials P„ may 
be easily determined and expressed in terms of f3j and e v upon taking (ITS]) into account. 
As an illustration we list here the first four polynomials 

4e r (l + z/)(2 + i/ 



AO 



7 + 2z/ +Pj 

4e r (16 + i/(13 + 3v))^- 
ll + 2z/ 



(21) 
(22) 



^(e r ,)9i) 

8e r (30 + z/(17 + 3z/))/?| 
15 + 2i/ 



48e r2 (l 



i/) (2 + i/) (5 + i/) (6 + i/) 



15 + 2i/) 



(23) 



Pi i} (^ r 



16e r2 (2944 + i/(3404 + i/(1451 + 5z/(50 + 3i/))))^- 



A') 



'19 + 2z/) 2 



+ 



40£ r (16 + z/(7 + z/))/?| 



(24) 



19 + 2z/ 

The values of e r , f3i, /3 2 , for which a closed-form solution of a given parity i/ exist, are 
determined by solving the algebraic equations Pi 1 \e r ,(3i) = and Pn\e r : (3 2 ) = 0, 
simultaneously with the imposition of (3\ + (3 2 = 2 at a chosen value of n. The solution 
may be explicitly determined by calculating the coefficients in the functions g(r]j, (3j) 
via (120!) for j = 1,2. We may notice that for a fixed confinement potential with a specific 
value of the frequency u x determined by (I18|) . there is only one solution that is exactly 
known. This is a consequence of the fact that the energy e r enters equation (fl8|) as a 
coefficient of the quadratic term in the potential, which is in difference with the case of 
the sextic oscillator Schrodinger equation, where several solutions may be analytically 
determined at specific values of the anharmonicity parameter [18] . 

We consider here the lowest exact solutions of for the 2D asymmetric QD generated 
through setting n = 1,2, 3, 4, etc. In the simplest case of v = and n = 1, we have 
the equations P{ (e r , (3j) = — ^-+/3? = 0, which upon condition (3i+(3 2 = 2, are solved 
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by Pj = 1, e r = 7/8, corresponding to the functions g^°'(r]j, 1) = e~32"(l + y) that are 
nodeless. This yields the normalized exact solution of (JE]) in cartesian coordinates in 
the form 

ip (m (x,y) = 1 e-g-£(l + ^x 2 + p), (25) 

2^(12^ + r[-i] 2 + 2r[i] 2 ) 4 

which corresponds to the ground-state. Further examples obtained in this way, for both 
the ground and excited states, are given in Tabled] with the corresponding values of u x , 
\S\, and e r . In Figures [2] and [3j we have plotted the normalized wave functions ip±, that 
correspond to the pairs of degenerate states, obtained for uj x = 1/32 and uo x = 1/16, 
respectively 



Table 1. Confinement frequencies, w x , and the corresponding exact values of the 
relative motion energy, e r , in the case of u> y : lu x = 2 at a given parity v for several 
values of n. 



n 


V 




spin 


symmetry parity 


1*1 


£ r 


(wij^a) 


1 





1 

8 


(s) 


(+,+) 





7 

8 


(0,0) 


2 




1 

61 


(s) 


(+>+) 





11 

64 


(0,0) 






1 

16 


(s),(t) 


(+,+),(+,-) 


1 


11 

16 


(2,0) 


3 




Bo(5 + 2x/5) 


(s) 


(+,+) 





|(5 + 2 VE) 


(2,2) 






i£o(5-2>/5) 


(s) 


(+,+) 





|(5 - 2V§) 


(0,0) 






240(5 + ^) 


00,00 


(+,+),(+,-) 


|(V5 + 1) 


^(5 + ^5) 


(0,4) 






240(5"^) 


(s),(t) 


(+,+),(+,-) 


KV5-1) 


^(5-V5) 


(2,0) 


1 


1 


1 

21 


(t) 


(->+) 





3 
8 


(1,1) 


2 




1 
128 


(t) 


(->+) 





13 
128 


(1,1) 






1 

32 


(s),(t) 


(--),(-+) 


1 


13 
32 


(3,1) 
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Figure 2. The exact wave functions tp(±) corresponding to the degenerate level of the 
relative motion energy e r = || 





Figure 3. The exact wave functions ip(±) corresponding to the degenerate level of the 
relative motion energy £ r = 

4. Frobenius method 

Now we develop a numerical scheme based on the Frobenius method (FM) for obtaining 
the approximate energies of an asymmetric QD ([8]) in the integrable case of u y : u x = 2. 
The method consists of expanding the solution of a differential equation into power series 
|19j . and was originally applied by Barakat and Rosner [18] to compute the spectrum 
of one-dimensional quartic oscillator confined by impenetrable walls at x = ±R. The 
energy eigenvalues of the system have been obtained numerically as zeros of a function, 
calculated from its power series representation. Moreover, it has been shown that the 
bound-state energies of the confined system approach rapidly those of the unconfined 
oscillator for increasing R [13], [HI US]- Here we show that solving the two separated 
equations f|T2|) with the help of the power series representation f|T5|) simultaneously with 
the condition f3\ + fa — 2 enables us to determine effectively the energy levels of the QD 
to very high precision. With the boundary conditions imposed on the functions gj{rjj) at 
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the points rjj = R, namely g(R,/3j) = 0, the above problem reduces to finding common 
zeros [5, e r ] of the functions g(R, 1 + 5) and g(R, 1 — 5). The functions g(R, 1 + 5) can be 
calculated numerically with an arbitrary accuracy from the power series representations 

K 

s(^) = E 6 FV l+I/ > ( 26 ) 

i=0 

truncated at a suitably high order K. The obtained approximations to g(rjj,(3j) are 
polynomials in the variables e T and 5, which can be taken to advantage in determining 
the numerical values of their common zeros. Since the exact bound-states of the QD 
correspond to the functions g(r]j,{3j) that vanish at 7] 3 • — >• oo, we expect the solution of 
the parabolically confined problem to approach the desired solution of the unconfined 
system, as the confinement range R increases. For the demonstration of convergence, 
we apply the above scheme to the QD with uj x = 1/16, for which the exact solution is 
known (see Tabled]). Table [2] shows how the numerical values [5, e£ M ], obtained for the 
confined system, approach the exact ones, as the truncation order K in the power series 
(]26p and the confinement range R increase. 

Table 2. Convergence of the energy eigenvalue and the corresponding separation 
constant determined by the FM. The underlined numbers agree with the exact results. 



K 


R 




t FM 


5 


60 


5 


0.6£ 


^74432520 


0.999811134 


70 




0.6£ 


^75024844 


1.000007959 


90 


5.5 


0.6£ 


^75000837 


1.000000296 


100 




0.6£ 


^75000006 


1.000000002 


120 




0.6^ 


^75000008 


1.000000003 


140 


6 


0.6£ 


^75000000 


1.000000000 



We have also applied the FM for determining the precise values of energies for the 
states that are not amenable to exact solution at various u y : u x — 2. The stability 
of numerical results was achieved by increasing K and R until the energy e v FM and 
the corresponding separation eigenvalue 5 stay fixed to the quoted accuracy, which 
determines the energy of the relative motion of a QD to that accuracy. The bound-state 
energies determined by the FM are listen in Table [3] up to 8-digit accuracy at particular 
values of ou x . The states for which the exact solution exists are marked bold, their 
energies can be shown to agree with the exact values to practically unlimited accuracy. 
The plot of ground state energy in function of Lnu;;,. is displayed in Fig. HI where we have 
marked a few of exact values of energy, which correspond to the closed-form solutions 
discussed in the previous Section. 

We may notice that the other integrable case, that of circular symmetry, where 
the problem reduces to the radial equation for the relative motion, can be also treated 
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Table 3. Low- lying bound-state energies at various confinement frequencies u x , as 
determined by the FM. The existing exact solution are marked bold. 





r(+,+) 


!"(->+) 


r(+, +)(+,-) 


r(-, -),(-,+) 


t FM 




L FM 


L FM 


1/64 


0.17187500 


0.172678376 


0.220586977 


0.225979324 


1/32 


0.293674143 


0.297716776 


0.386949414 


0.40625000 


1/24 


0.367590598 


0.37500000 


0.490314824 


0.521041536 


1/16 


0.505362736 


0.521827040 


0.68750000 


0.743921450 


1/8 


0.87500000 


0.931629324 


1.241576361 


1.386478626 


1/6 


1.100931183 


1.191668135 


1.595014888 


1.803819570 


1/2 


2.681851499 


3.142358616 


4.266819742 


5.035425098 


1 


4.773351606 


5.921649105 


8.099126926 


9.762556072 


2 


8.623556558 


11.317113654 


15.569391839 


19.083841664 



within the FM. In that case the solutions has a generalized power series representation 
and the highly accurate values of energies can be determined with a small computational 
effort by the scheme developed recently in our previous work [T6] . 



1.2 
1 

. 8 
% 0.6 
. 4 
0.2 



-8 -7 -6 -5 -4 -3 -2 
Ln [w x ] 

Figure 4. The relative motion ground state energy e r of a QD as a function of 
hiaw x (ujy : u> x = 2) determined by the FM. A few values of energy, for which closed- 
form solutions exist are marked with crosses. 

On the example of a QD with confinement frequency = ^, we use the exact 
ground-state energy and the very precise numerical values for higher bound-states 
determined by the FM (see Table |3]) to check the convergence of the popular numerical 
methods applied directly to the two dimensional Schrodinger equation (jHJ). First, we 
consider the Rayleigh-Ritz (RR) method with the basis set of the two dimensional 
harmonic oscillator eigenstates, calculating the approximate values of energies by 




Quasi- exact solutions for two interacting electrons 



12 



Table 4. Convergence of low-lying states energy eigenvalues calculated by the 
Raylcigh-Ritz method for a QD with lo x = 1/64. The underlined numbers agree 
with the exact result or that determined by the FM. 



D 



£ 



r(+,+) 
RR 



5 



*(-,+) 
RR 



£ 



r(+,+) 

RR 



£ 



r(+,-) 

RR 



■RR 



£ 



r(-,+) 

RR 



25 0.1718 911 0.1726 898 0.220 6744 0.220 6004 0.2259 834 0.2260368 

64 0.1718 856 0.1726 830 0.220 6472 0.2205 919 0.2259 803 0.2259 999 

144 0.1718 823 0.1726 804 0.220 6274 0.22058 91 0.225979 6 0.2259 878 

256 0.1718805 0.1726795 0.2206169 0.2205881 0.2259795 0.2259839 



Table 5. Same as in Table 0] but calculated by the DVR method. 



e DVR 



256 


30 


0.1719120 


0.1727763 


0.2206480 


0.2214837 


0.2260438 


0.2277602 


400 




0.1719294 


0.1727757 


0.2206258 


0.2215710 


0.2260337 


0.2277544 




35 


0.1718335 


0.1726833 


0.2204264 


0.2206146 


0.2259909 


0.2260835 


576 




0.1718460 


0.1726824 


0.2204845 


0.2205998 


0.2259843 


0.2260752 


900 




0.1718557 


0.1726815 


0.2205324 


0.2205924 


0.2259816 


0.2260703 




40 


0.1718483 


0.1726799 


0.2204538 


0.2205956 


0.2259819 


0.2259880 



separate diagonalization of the matrix Hamiltonian corresponding to the eigenfunctions 
with definite (x, y)— parity in a manner similar to that used in Ref. [11]. In Table H] 
the obtained values of energies e RR are presented for increasing dimension D of the 
RR matrix. The same QD example is used for testing the DVR (discrete variable 
representation) method with the trigonometric basis for the square grid on a square 
whose sides have length R; Xi = y { = —R + 2Ri/N, i = 1, N — 1 (for details see [20]). 
The numerical values of energies s^vr obtained by diagonalization of the Hamiltonian 
matrix in point grid representations with dimension d = (N — l) 2 , are presented in Table 
[5] for several values of R and d. We may see that in both considered methods, the rough 
values of energies can be determined with the use of small basis sets but the convergence 
is rather slow although the computational effort increases quadratically with the size of 
the basis. 

5. Concluding remarks 

We discussed the case of a QD consisting of two interacting electrons confined in an 
anisotropic harmonic potential V(x,y) = m^iulx 2 + ojyy 2 )/2. In the case of uj v : cu x = 2 
we developed a numerical scheme based on the Frobenius method that allows an efficient 
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determination of the spectrum. Moreover, we have shown that a closed-form analytical 
solutions for the relative motion can be obtained at particular values of uj x . We 
illustrated using both the exact and the very precise values of energies determined by FM 
for testing the convergence of the Rayleigh-Ritz and DVR approximation methods. We 
have shown that both the methods determine well the approximate values of low-lying 
states energies but their convergence rates are quite slow. 

The exact values of energies and the corresponding wave functions obtained in our 
work can be useful for investigating the performance of other approximation methods 
in the case of the asymmetric confinement. They can be also used to build the exact 
density functional and testing the possible ways of constructing approximations to it. 
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